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Abstract 

We investigate Wilson twisted mass fermions in the quenched approxima¬ 
tion using different definitions of the critical bare quark mass rric to realize 
maximal twist and, correspondingly, automatic 0(o) improvement for phys¬ 
ical observables. A particular definition of rric is given by extrapolating the 
value of rric obtained from the PGAC relation at non-vanishing bare twisted 
quark mass /x to // = 0. Employing this improved definition of the critical 
mass the Wilson twisted mass formulation provides the possibility to perform 
reliable simulations down to very small quark masses with correspondingly 
small pion masses of ~ 250 MeV, while keeping the cutoff effects of O(a^) 
under control. 


1 


1 Introduction 


The combination of an automatic 0(a)-improvement, the infrared regulation of 
small eigenvalues and fast dynamical simulations render the so-called twisted mass 
fermions [1112113 a most promising formulation of lattice QCD. There are already 
a number of results for the quenched approximation in 0 El- Also full QCD sim¬ 
ulations with this approach have been performed and revealed a surprising phase 
structure of lattice QCD [aiH]. On the theoretical side, extensive calculations in 
chiral perturbation theory were done jnmniiii]. 

One of the problems that appeared in the discussion of twisted mass fermions is 
the exact dehnition of the critical bare quark mass. In jT2] it was shown that using 
the critical value of the bare quark mass where the pion mass vanishes, leads to a 
quark mass dependence of the pion decay constant /^r that deviates strongly from 
the expected linear behavior. This ’’bending phenomenon” was observed when the 
quark mass rriq obeys the inequality 

arriq < . ( 1 ) 

The bending phenomenon was interpreted in a way that for such small values of the 
quark mass the Wilson term becomes the dominant term in the lattice Wilson-Dirac 
operator, thus leading to large lattice artefacts. 

In refs. [mill it has been suggested to take the value of the critical mass where 
the PCAC quark mass vanishes. In chiral perturbation theory it was demonstrated 
that this dehnition leads to a stable situation of ’’maximal twist” and hence to a 
reduction of lattice artefacts that otherwise would be enhanced at very small values 
of the quark mass. In refs. [m nni the theoretical background of such a dehnition 
is discussed on a general ground. 

In this paper, we employ the dehnition of the critical mass by computing at non¬ 
zero twisted mass the critical hopping parameter as a functions of the twisted mass 
parameter and hnally taking the limit to zero twisted mass. Such a dehnition is ex¬ 
pected to eliminate large O(a^) artifacts at small values of the pion mass and should 
hence avoid the bending phenomenon. In the following we will shortly describe the 
encouraging results employing this improved dehnition of the critical mass. A more 
detailed scaling analysis of mesonic quantities will be reported elsewhere HZI. 


2 Wilson twisted mass fermions 

In this paper we will work with Wilson twisted mass fermions that can be arranged 
to be 0{a) improved without employing specihc improvement terms 0 |2]- The 
Wilson tmQCD action in the twisted basis can be written as 

S'[f/,'0,'0] = + mo + , (2) 
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where the Wilson-Dirac operator Z)w is given by 


^ + V/,) - aV* V^] 

Ai=0 ^ 


( 3 ) 


and and V* denote the usual forward and backward derivatives. We refer to 
1121 for further unexplained notations. The dehnition of the critical mass rric will be 
discussed in detail in the next section. 

We extract pseudoscalar and vector meson masses from the correlation functions 
at full twist (mo = rric): 

C“(xo) = a = 1,2 

X 

CHxo) = 0 = 1,2 

k=l X 

Ct(xo) = 

k=l X 


where we consider the usual local bilinears = -075^'0; — '^lklb'^4’ and 

= i’aoki'^. 

The untwisted PCAC quark mass mpcAc can be extracted from the ratio 


W^PCAC — 


E,(ao2iS(x) ^“(0)) 


o = 1,2 . 


(4) 


Using the PCAC relation we can also compute the pion decay constant at max¬ 
imal twist without requiring any renormalization constant (see naa): 

A=(^l(0|P»| a =1,2. (6) 


3 Definition of the critical mass 

The Wilson tmQCD action of eq. (j21) can be studied in the full parameter space 
A special case arises, however, when mo is tuned towards a critical bare 
quark mass me. In such, and only in such a situation, all physical quantities are, or 
can easily be, 0{a) improved. The critical quark mass, or alternatively the critical 
hopping parameter Kc = 2am +8 ’ hxed in the actual simulation to 

achieve automatic 0(a)-improvement at maximal twist. 

In general, the dehnition of the critical mass has an intrinsic uncertainty that 
comes at 0{a) for Wilson fermions. This could be improved to an O(a^) uncer¬ 
tainty by using Clover fermions. If a dehnition of Kc from a linear extrapolation of 
{rriT^a)'^ —> 0 is used, one has in addition an unrelated systematic error coming from 
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the long way of extrapolation from large pion masses of 0(600 MeV). 

A better definition of Kc can be obtained by the following procedure: 

• At hxed non-zero twisted mass parameter the hopping parameter Kc{fJ^a) is 
determined as the point where the PCAC quark mass of eq. m vanishes. The 
non-zero value of the twisted mass allows a safe interpolation in this case. 

• As a further step, an only short linear extrapolation of Kc{fJ^a) from small 
values of pa to pa = 0 yields a dehnition of Kc which is expected [THl Uni to 
lead to small lattice artefacts, in particular also at small quark masses much 
below values indicated in the inequality of eq. (P). 

In hgure ^ we show an example of this extrapolation at (3 = 5.7, where we also 
indicate the Kc-value determined by the vanishing pion mass squared. 



Figure 1: Determination of the critical mass: I/kc versus pa at /? = 5.70, ex¬ 
trapolation to pa = 0, the open diamond indicates the 1/kc value determined by 
{rriT^ay —> 0 at p = 0 for unimproved Wilson fermions. 

The possible choices for the determination of the critical bare quark mass were 
also discussed in chiral perturbation theory [nnn. Recently, a dehnition of maxi¬ 
mal twist from parity conservation has also been investigated in jS]. There, however, 
the critical masses mc(pa) were not extrapolated to pa = 0, but were used at the 
respective twisted mass parameter at which they were determined. 
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4 Numerical results 


In this section we will provide a comparison of Wilson twisted mass results for the 
pion mass, the pion decay constant and the p-mass obtained with two different 
definitions of Hc- The hrst, “pion mass” definition, refers to the limit {rn-,^aY —> 0 
taken at pa = 0. The second, “PCAC quark mass” dehnition, corresponds to the 
determination of Kc(pa) at fixed fia ^ 0 from the vanishing of the PCAC quark 
mass, ampcAcif^o,) and then taking the limit > 0). The latter Kc values 

were determined on a subset of 0(150) configurations for three couplings, namely 
d = 5.70, 5.85 and 6.0. The simulations were performed for a number of bare quark 
masses in a corresponding pion mass range of 250 MeV < m-,^ < 1200 MeV using 
a multiple mass solver BSI on 0(400) gauge field configurations generated with the 
Wilson plaquette gauge action. 

In general the time component of the tensor correlator provides a better signal 
for the vector meson mass compared to the axial-vector correlator at maximal twist. 
This observation was made independently in [Oj and is manifest in the smaller error 
bars displayed in figures of this section. 



Figure 2: Comparison of results at /3 = 6.0 for the pion mass for different definitions 
of the critical mass. 

The first quantity we investigated was the pion mass. In figure |21 we show the 
behavior of {m.^aY versus pa for the two definitions of Kc- The results for the PCAC 
quark mass definition show a linear behavior down to very small bare quark masses. 
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while a small non-zero value is approached in the chiral limit with the pion mass 
definition of Kc- This residual 0{a) pion mass at /io = 0 can be attributed to the 
0{a) error in the critical mass determined by the vanishing of the pion mass in the 
pure Wilson case. 

At = 5.7 (lattice size 12^ x 32, a ~ 0.17 fm) we computed correlation functions 
also for a third, intermediate definition of the critical mass: we chose the point where 
the PCAC quark mass vanishes at small, but non-zero value of the twisted mass. 
In figure El we show the results for the pion decay constant and the vector meson 
mass. It is obvious that already this intermediate definition of Kc substantially 
diminishes the bending tendency of the data when approaching the chiral limit. Such 
a definition of the critical hopping parameter fi:c(/ro « 7 ^ 0) seems to be sufficient to 
eliminate the deviation from a linear behavior for all twisted mass values /xa > /io o- 
The results for Kc{^ia = 0) from the PCAC quark mass definition show a straight 
behavior for both observables down to pion masses of about 250 MeV. The bending 
phenomena reported in IQ can thus be attributed to the pion mass definition of 
the critical mass in this simulation. Note that at /? = 5.7 we are working at a rather 
coarse value of the lattice spacing. 

At /3 = 5.85 (lattice size 16^ x 32, a ~ 0.12 fm) we can now compare our previous 
and new Wilson twisted mass results with the data from overlap fermion simulations 
IQ which were obtained on a smaller lattice volume (12^ x 24). Using the definition 
of Kc from the vanishing of the PCAC quark mass, the bending near the chiral limit 
vanishes almost completely both for the pion decay constant and the vector meson 
mass extracted from the tensor correlator (see figure |3)). In both cases the overlap 
fermion results lie slightly higher than the twisted mass fermion results. This fact 
might be explained by residual finite volume effects and/or different lattice artefacts. 

Finally we show in figure El the results aX j3 = 6.0 (lattice size 16^ x 32, a ~ 
0.09 fm). Although the physical volume is smaller than the one aX j3 = 5.85 we are 
able to confirm the absence of bending for small masses when using the PCAC quark 
mass definition of Kc- For the pion decay constant we observe a very good agreement 
with the overlap fermion results of ref. 120] and with the twisted mass results using 
the parity conservation definition of the critical mass from ref. lUj. For the vector 
meson mass we refrain from comparing with the parity conservation definition, since 
the error bars would cover both our results with the pion mass and PCAC quark 
mass definition of the critical mass. 


5 Conclusions 

We have studied Wilson twisted mass fermions with different definitions of the crit¬ 
ical mass. The ’’bending phenomena” in basic observables like the vector meson 
mass and the pion decay constant using the critical mass from vanishing pion mass 
at /m = 0 reported in IQ is absent for the choice of Kc from the vanishing of the 
PCAC quark mass in the limit /ra —*■ 0. This is true even at small values of the 
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pion mass close to the experimentally measured value where maybe contact to chiral 
perturbation theory can be made. 

In ref. jH], the same dehnition of the critical mass has been employed, but at 
fixed twisted mass parameter separately at each simulation point. The results of 
this reference taken together with the hndings in the present paper lead to the 
conclusion that the PCAC quark mass dehnition of the critical hopping parameter is 
a crucial element of twisted mass simulations in order to keep O(a^) effects well under 
control at small quark masses. This conclusion is in accordance with theoretical 
considerations [Ml 13 ESI EE!. 

Employing the PCAC quark mass dehnition of the critical mass the Wilson 
twisted mass formulation provides the possibility to perform reliable simulations at 
very small quark masses ~ 250 MeV). Thus, Wilson twisted mass fermions can 
be used to explore really light quarks on the lattice - as light as with overlap fermions, 
but at considerably lower cost (see |2I] for a detailed cost comparison). This opens 
a very promising prospect for dynamical fermion simulations and renders twisted 
mass fermions as a real alternative to staggered fermions without the conceptual 
difficulties of the latter. 
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Figure 3: Comparison of results at /3 = 5.70 for the pion decay constant (top) and 
vector meson mass (bottom) for different definitions of the critical mass. 
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Figure 4: Comparison of overlap fermion and Wilson twisted mass results aX (3 = 5.85 
for the pion decay constant (top) and vector meson mass (bottom) for different 
definitions of the critical mass. 
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Figure 5: Comparison of results at /3 = 6.0 for the pion decay constant (top) and 
vector meson mass (bottom) for different definitions of the critical mass. 
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